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HEAT TRANSFER IN THE TURBULENT INCOMPRESSIBLE BOUNDARY LAYER 
IV - EFFECT OF LOCATION OF TRANSITION AND PREDICTION OF 
HEAT TRANSFER IN A KNOWN TRANSITION REGION 
By W. C. Reynolds, W. M. Kays, and S. J. Kline 


SUMMARY 

The effect of the location of transition on the heat transfer to the 
turbulent incompressible boundary layer is analyzed. The analysis indi- 
cates that considerably higher heat-transfer rates may occur for some 
distance downstream if the transition is very late. The results of a 
limited experimental investigation are in substantial agreement with the 
results of the analysis. 

If the extent of the transition region is known, the analysis also 
allows adequate prediction of heat-transfer coefficients within the tran- 
sition region. The nature of this analysis is such that it should pre- 
dict local shear coefficients in the transition region equally well. 


INTRODUCTION 

This report is the last of a series of four covering a three-year 
investigation of heat transfer through the turbulent incompressible 
boundary layer with arbitrary wall-temperature variation (ref. l) . The 
first report describes the experimental apparatus and presents results 
of experiments with constant wall temperature (ref. 2) . The second re- 
port gives the results of experiments and analyses for a step temperature 
distribution (ref. 3). In the third report the step-function analysis is 
used as the basis for predicting heat-transfer rates in several cases 
where the plate is nonisothermal (ref. 4) . These predictions are com- 
pared with experiments, and good agreement is founds and a simple method 
for handling problems where the wall temperature varies in an arbitrary 
manner is presented. In the present report an analysis of the effect of 
the location of transition on heat transfer in the turbulent boundary 
layer downstream is presented and compared with experiments. 

Experiments have shown that delayed transition from a laminar to a 
turbulent boundary layer results in higher heat-transfer rates in the 
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turbulent boundary layer than would be obtained if the boundary layer 
were turbulent from the leading edge (e.g., ref. 5). This is believed 
to be due to two effects. The first is a hydraulic effect; that is, the 
turbulent boundary layer behaves as if it had originated at some point 
downstream of the leading edge of the flat plate. The second is a 
thermal effect; that is, the heat-transfer rates in the laminar layer 
preceding the turbulent layer are much lower than those that would occur 
over the leading portion of the plate in a turbulent layer. Thus, the 
energy in the turbulent layer after transition is considerably less than 
it would have been if the layer had been turbulent from the start. Be- 
cause of the combination of these effects, considerably higher heat- 
transfer rates may occur in the turbulent part of the boundary layer for 
some distance downstream if the transition Reynolds number is high. 

An attempt has been made to predict analytically the effect of a 
late transition, based on the notion of separate contributions of hydrau- 
lic and thermal effects. The limited amount of experimental data ob- 
tained in this investigation for a variable transition Reynolds n umb er 
is correlated quite adequately on the basis of this analysis. However, 
it was not possible to delay transition to Reynolds numbers above 400,000 
in the available wind tunnel. Consequently, complete confirmation of the 
analysis must await the availability of experimental data for higher tran- 
sition Reynolds numbers. 

Since the length of the transition region is generally about the 
same as the length of the laminar region, estimations of heat transfer in 
the transition region are often desired. If the extent of the transition 
region is known, the local Stanton numbers in the transition and turbu- 
lent regions may be determined by employing the abrupt-transition anal- 
ysis as a point concept. This concept assumes that the flow at a given 
point is either laminar or turbulent and that the average condition of 
the flow at a prescribed distance from the leading edge can be represented 
by the value of an intermittency factor that prescribes the fraction of 
the flow that is turbulent. The known expressions for laminar and tur- 
bulent heat-transfer coefficients are then assumed to hold at each point, 
and the mean coefficient is found as a weighted average. The heat- 
transfer coefficients found in this manner agree well with the data ob- 
tained with natural transition in the present investigation. It is be- 
lieved that the method should work equally well for predicting local shear 
coefficients in the transition and turbulent regions. 

This investigation was carried out at Stanford University under the 
sponsorship and with the financial assistance of the National Advisory 
Committee for Aeronautics. 
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SYMBOLS 

location of start of transition, ft 
location of end of transition, ft 
local friction factor, T w /(pu^/2) 
laminar friction factor 
turbulent friction factor 

specific heat at constant pressure, Btu/ (lb) (°F) 

St turb( Re x> Re sV st turb( Re x>‘°) 
free-stream mass velocity, pu^, lb/ (hr) (ft) 

local convective heat-transfer coefficient, q"/At, 

Btu/ (hr) (sq ft) (°F) 

thermal conductivity of fluid, Btu/ (hr) (ft) (°F) 
unheated starting length, ft 

probability distribution function for location of transi- 
tion, P(Re g ) 

Prandtl number, pCp/k 

heat flux at wall, Btu/ (hr) (sq ft) 

Reynolds number at abrupt transition 

Reynolds number based on x, Gx/p 

Reynolds number based on co 

location of virtual origin, ft 

local Stanton number, hj/GCp 

laminar Stanton number 

local Stanton numbers on isothermal plate for turbulent 
boundary layer after abrupt transition at s 


turbulent Stanton number 


location of abrupt transition, ft 
average location of abrupt transition, ft 
absolute wall temperature , °R 
absolute free-stream temperature, °R 
t - t , °P 

At for fictitious part of turbulent layer, °F 
mean temperature of heated strip, °F 

tm " too, °F 

wall temperature, °F 

free-stream temperature, °F 

velocity in boundary layer, ft/ sec 

free-stream velocity, ft/sec 

flow rate, Ib/hr 

distance from leading edge, ft 

distance from wall, ft 

coordinate normal to x and y, ft 


fraction of boundary layer at x that is turbulent 
P(Re x ) 


boundary-layer thickness , ft 

laminar -boundary-layer thickness at a, ft 

turbulent -boundary-layer thickness at b, ft 


momentum thickness* r foo/ujll - (uKi/Oldy, f 



5 


(i viscosity, lb/ (hr) (ft) 

| distance from virtual origin, ft 

p fluid density, lb/cu ft 

a variable of integration 

x shear stress at wall, lb/sq ft 

w 

co standard deviation of s from mean transition-point loca- 

tion, ft 


TRANSITION MODEL 

Although a great deal is known about the heat-transfer and skin- 
friction characteristics of laminar and turbulent boundary layers, very 
little is known about heat transfer in the region of transition that ex- 
ists between the turbulent and laminar regimes, or about the effect that 
the location of the transition region has on heat transfer and skin fric- 
tion in the turbulent layer after transition. In order to estimate the 
effect of transition-point location on turbulent heat transfer for en- 
gineering design calculations, one of the following three assumptions 
has usually been made: (l) The turbulent layer behaves as if it had 

started at the leading edge; (2) it behaves as if it had started at a 
"transition point;" or (3) the turbulent layer after transition has the 
same thickness as the laminar layer before transition. The heat-transfer 
rates and friction factors determined by these approximate treatments 
will not be too much in error, if the point under consideration is very 
far downstream of the transition region and if transition ends at rela- 
tively low Reynolds numbers. In some applications, however, the transi- 
tion may be quite late; not infrequently it extends to Reynolds numbers 
in excess of 10 , and in such a case the late transition may cause a con- 
siderable increase in the turbulent skin friction and heat transfer for 
a surprisingly large distance downstream of transition. A realistic 
analysis of the transition effect is therefore needed; the analysis should 
be supported by experimental data, since assumptions will be required by 
the nature of the problem. 

In making an analysis of the transition effect, it is important that 
a reasonable model of transition be employed. The recent experiments of 
Schubauer and Klebanoff (ref. 6) indicate that the laminar boundary layer 
on a flat plate becomes turbulent quite suddenly at various "spots" on 
the plate. The location of these spots of "abrupt" transition varies with 
time, but the amount of the boundary layer that is turbulent at any dis- 
tance downstream of the leading edge seems to depend only on the distance. 
These observations suggest a model of the transition phenomena that should 
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be suitable as the basis for an analysis of the effect of the location 
of transition on heat transfer and skin friction. This model is shown 
by figure 1. The flow is divided into "slabs/' and the flow in each 
"slab" is assumed to be independent of the flow in the other slabs. In 
each slab, transition occurs abruptly at some point, and the location of 
the "local abrupt transition" differs from slab to slab. However, if 
the slabs are made very thin, the locus of the transition points becomes 
a continuous curve, as is indicated by figure 1. It will be assumed that 
this curve is invariant in time. This assumption is in contradiction 
with the observations of Schubauer and KLebanoff, but since a statistical 
approach will be taken in the analysis no serious error should result. 

Well downstream of the transition locus, the effect of the variation in 
the location of the transition point will be smoothed out, and the local 
Stanton numbers and friction factors should approach the values for a 
boundary layer that is turbulent from the leading edge. 

By examining this model statistically, it is possible to derive ex- 
pressions for the local heat-transfer rate and skin friction that apply 
over the entire length of the plate. The Stanton number may be expressed 
in terms of the Stanton number for a laminar boundary layer, the Stanton 
number for a turbulent boundary layer undergoing an abrupt two-dimensional 
transition, and the fraction of the flow that is turbulent at any point 
x downstream of the leading edge. A similar expression can be derived 
for the local friction factor. 

The location of transition in any slab will be denoted by s; s 
may be described by a probability distribution P(s) where P(s)ds 
represents the probability that transition occurs between s and s + ds. 
The probability that, at any point x, the flow in a given slab will be 
turbulent will be denoted by 3 (x) , the intermittency factor; 3 (x) is 
also the fraction of the boundary layer that is turbulent at x. The 
fraction of the boundary layer that is laminar at x is taken as 
1 - 3(x). Since S(x) is the probability that transition has occurred 
before x, 3(x) is related to P(s) in the following manner: 



Moreover, since the boundary layer eventually is entirely turbulent, 

( 2 ) 

The mean location of transition is defined by 



3(oo) = 1 = f p(s)ds 


( 5 ) 
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where s m represents the average location of transition. Note that if 
P(s) is symmetric, then s m will be the symmetry point; and, for approx- 
imately normal distributions, the maximum value of P(s) will be P(s m ). 

The rate of heat transfer to a turbulent slab of the boundary layer 
will depend on the usual flow parameters, the location of the point in 
question (x) and the location of the abrupt transition (s). This depend- 
ence may be indicated by writing 

Iturb = ‘Iturb^ 3 ) ( 4 ) 

The average heat-transfer rate for all of the turbulent slabs at any lo- 
cation x is defined as Since q" (x;s)P(s)ds is the contri- 
bution to the heat flux at location x due to all slabs that went tur- 

bulent between s and s + ds, it follows that the total contribution 
to the heat flux at x from all slabs that are turbulent is 



Similarly, the average heat-transfer rate for all of the laminar slabs 
at any location x is 

^Zam = q Zam^ x ^ 


The total local heat-transfer rate at x for both the laminar and tur- 
bulent portions of the flow is therefore the sum of the laminar and tur- 
bulent contributions: 

q;(x) = [l - P to] 1^00 + P(x)i^ b (x;s) (?) 

By dividing by suitable constants and using equations (l) and (5) , this 
relation may be written as 


-Re x 

St(x) = [l - p(x)]st lam + £ St turb (Re x ;Re s )P(Re s )dRe 


( 8 ) 


Similarly, for the friction factor, 

.Re. 


C f (x) = [l - P(x)] C fflam + C f ^ turb (Re x ;Re s )p(Re s )dRe £ 


( 9 ) 
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Therefore, if the local Stanton numbers for a turbulent boundary layer 
after an abrupt transition are known, and if the fraction of the flow 
that is turbulent at any point x is known, the local skin friction an d 
heat transfer for the nonabrupt transition can be determined by perform- 
ing the integrations of equations (8) and (9). It should be noted that 
this statistical approach allows calculation of the local Stanton numbers 
and friction factors in all the laminar, transition, and turbulent regions 
of the boundary layer. 


p(x) for "Abrupt" Transition 

If the boundary layer is tripped at a relatively low Reynolds num- 
ber in such a way that transition occurs abruptly at the same value of 
x across the entire plate, p(x) is simply a "jump function," and 

p(x) = 0 x < s 

p(x) =1 x > s 

In this case the probability density is discontinuous (a "pulse" of unit 
area), and the integrals must be evaluated in the Stilt jes sense, which 
gives 


St(x) = St Zam (Re x ) 

x < s 

St(x) = St turb (Re x ;Re g ) 

x > s 

C f (x) = C f)lam (-Re x ) 

X < s 

Cf( x ) = C'-p^f 11 - r -h(Re- x ;Re g ) 

X > s 


where s is the location of the abrupt transition. 


p(x) for "Natural" Transition 

The data of Schubauer and KLebanoff (ref. 6) indicate that, for 
natural transition, p resembles an error integral. If the parameters 
of the error curve were known, the integrals could, in principle, be 
evaluated. Unfortunately, the form of the abrupt-transition Stanton num- 
ber expression is such that the integrations may not be performed in 
closed form. However, for values of x > s m the integrals may be approx- 
imated by replacing s by s m in the expression for the abrupt- 

transition Stanton number. The justification of this is demonstrated by 
the following sketch: 


CX-2 
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St(Re : 



-St(Re x ;0) 


This sketch shows the values of the Stanton numbers that would obtain 
at x for various locations of the abrupt transition. If x exceeds 
s m , then St(Re x ;Re s ^) lies between St(Re x ;0) and St(Re x ;Re x ) , and a 

reasonable approximation is 



R ^turb ^ Re x > Re s ) p ^ Re s ^ *® e s 


r R e x 

J st turb( Re x^ Re s m ) p ( Re s ) dRe , 

0 m 

/• Re x 

^°turb v,i ' c x’ ric s m ‘' J P\,Re g )dRe 

m 0 

p(Re x )St turb (Re x ;Re ) x > s. 


m 

(10a) 


For values of x less than s m this approximation has no meaning, and 
thus another approximation must be used. For x < s m , the quantity 
Stturb( Re x> Re s) be approximated by 

St(Re x ;Re g -aB ) 
m 


where co is the standard deviation of s from the mean transition-point 
location. This approximation will be used in the range (s m - 2oo) < x < s m , 

and leads to 



St turb^ e x^ Re s) P C Re s) dRe s 


3(x)St tur b( Re x^ Re s m - Re 2cu) 


(s m - 2cn) < x < s m (10b) 


10 


In the laminar region (i.e., for x < (s m - 2oo) ) , (B is small (p < 0.03), 

and thus the turbulent term may he neglected entirely. To summarize, 
then, the local Stanton number may therefore be approximated by 

St(x) = St Zam (Re x ) x < (s m - 2co) (lla) 

St(x) - (1 - P)St !m (Re x ) + SSt turl) (Re x J R. % - Re^,) (s n - Sol) < x < s m (111.) 
. (1 - P)St lsm (Se x ) + ’ s » (llc > 

Similar expressions for the local friction factor may be obtained: 

C fM = C f,Iam (Re x) x < ( s m " 2®) ( 12a ) 

C f (x) = (1 - P)C f) j am (Re x ) + PC f>turt (Re x ;Re Sm - Re^) (s m - 2co) < x < s m (12b) 

C f (x) = (1 - P)C f?Zam (Re x ) + pC f;turb (Re x ;Re Sm ) x > s m (12c) 

Assuming that the function p(x) is known, it is necessary only to 
obtain expressions for the local Stanton number and friction factor for 
a turbulent layer that underwent transition abruptly at Re s . An anal- 
ysis leading to such expressions is presented herein. 


ABRUPT-TRANSITION ANALYSIS 

In the abrupt-transition analysis it is assumed that the turbulent 
boundary layer behaves hydrodynamic ally as if it "started" at some "vir- 
tual origin" (see fig. 2(a)). The location of this virtual origin, which 
will be denoted by r, is determined by applying the momentum and conti- 
nuity theorems to the transition region (fig. 2(a)). The friction fac- 
tors for the turbulent layer may then be calculated from equations for 
boundary layers that are turbulent from the leading edge, if the Reynolds 
number is based on the distance x - r. The local Stanton numbers for 
the turbulent layer could be calculated in the same manner, except that, 
if the plate is isothermal, the total heat transferred to the laminar 
layer before transition exceeds what would be transferred to the ficti- 
tious turbulent layer upstream of transition, and consequently the con- 
servation of energy condition could not be satisfied for the transition 
control volume. However, by allowing the plate temperature to be higher 
in the fictitious portion of the turbulent layer, the energy condition 
can be satisfied. The "equivalent thermal problem" is shown in figure 
2(b). The local Stanton numbers for this equivalent problem may be deter- 
mined by employing nonisothermal heat-transfer theory (see refs. 3 and 4). 
This approach allows momentum, continuity, and energy conditions to be 
satisfied at the abrupt transition, and thus the analysis has a firm 
physical basis. 
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Equivalent Hydrodynamic Problem 

Consider a laminar boundary layer flowing over a flat plate as shown 
by figure 2(a). A transition region extends over a small distance dx 
from a, where the flow is purely laminar, to b, where a fully developed 
turbulent boundary layer exists. The laminar-boundary-layer thickness at 
point a will be denoted by and the turbulent -boundary-layer thick- 
ness at b by The local wall shear stress in the transition region 

is t^(x) . The turbulent boundary layer is assumed to behave hy dr ©dynam- 
ically as if it had originated at the virtual origin r; dx is small so 
that the shear stress may be considered constant. 


Application of the momentum theorem to the control volume aa'b'b 
yields 



+ u w f , , 
oo a ! b’ 



dy + dx 


(13) 


* 




Here w^-^, is "the mass-flow rate across the surface a'b' that enters 
the control volume with a velocity component in the x-direction of u^. 
Use of the equation of continuity for the same control volume gives 




/*" a r** 

/ P u b Z + UooP(&b - S a> + w a'b ' = / P u b l d y ( 14 ) 

°(aa') 0 (bb ' ) 

Multiplying equation (l4) by u m and combining with (13) give 

.&5 

'0 *" ~ •'O 


f*b f* a 

/ ( pu bZ u « " pu bZ )dy ~ J ( pu bZ U oo - P^bl)^ = T w ( 15 ) 


Using the definition of the momentum thickness 6 and dividing by pu 
reduce this to 


- 0 a - “*5 


a 2 

n 

~ 00 


(16) 


where 0 a and 6^ are the laminar and turbulent momentum thicknesses at 

a and b, respectively. Since dx is small, the term involving the 
shear stress is of higher order and may be neglected. Thus, equation 
(16) reduces to 




6 


a 


®b 


(17) 
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Equation (17) is the criterion used to define the virtual origin of the 
turbulent boundary layer. Both continuity and momentum are satisfied by 
the condition of equation (l7). 

The velocity profile in the laminar boundary layer may be approx- 
imated by (ref. 7 , p. 69) 

if (-> 

Thus, the laminar momentum thickness at point a may be evaluated as 
follows : 


e = s 

a a 


r (i ° - i ° 5 - i ° 2 + i ° 4 - i ° 6 ) ^ 


280 °a 


For the turbulent boundary layer in the range of Reynolds numbers found 
at transition, a good approximation to the velocity profile is 

(?.a\ 




Evaluation of the turbulent momentum thickness at b gives 


*> - *» 7 1 (■ 


C, 1 / 7 - a 2 / 7 )do 




Substituting equations (19) and (2l) into (17) yields 


^ = 1.433 
8 a 


The momentum equation of the turbulent incompressible flat-plate 
boundary layer may be -written as 


* /' (?)(■ - ?) 


Substitution of the laminar velocity profile (eq. (18)) into equation 
(23) leads to (ref. 7 , p. 69) 



a 


4.64Re" 0,5 

a 


(24) 
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For turbulent flow in the range 10 < Re x < 10', the local friction fac- 

tor can be given by the Blasius formula (see ref. 7, p. 117), as 



0.0228Reg 0,25 


(25) 


Substituting equations (20) and (25) into (23) gives an expression for 
the turbulent-boundary-layer thickness: 


| = 0.0376Re| 0 ' 2 


(26) 


Here | is measured from the virtual origin of the turbulent layer, 
point r.. Assuming that the turbulent layer after transition behaves as 
if it had originated at the virtual origin r, equation (26) then gives 


_&b_ 

b - r 


-O 9 

0.37GRe, u 
o-r 


(27) 


By letting dx become very small, a approaches b, and, in the limit, 
combination of (22), (24), and (27) results in 

Re r = Re s ^l - 36.3Reg 0,375 ^ (28) 

where Re s = Re a = Re-^. Equation (28) defines the location of the vir- 
tual origin of the turbulent boundary layer. 


Effect of Transition-Point Location on Local 
Turbulent Friction Factors 

As shown by reference 7 (p. 117), combination of equations (20), 
(23), and (25) results in the following relation for the turbulent fric- 
tion factor: 


C f (Re|jO) 

2 


0.0296Re|°* 2 


(29) 


14 


where again £ is measured from the virtual origin of the turbulent 
layer. Since the boundary layer behaves as if it begins at r, the 
friction factor at some point x may be determined from (28) and (29) 
as 


C f (Re x ;Re s ) 


= 0.0296Re 


- 0.2 

x-r 


- 0.0296Re x 0,2 1 - (l - 36. 3Reg 0 ' 375 ) "°‘ 2 


(30) 


The term before the bracket in equation (30) is the local friction coef- 
ficient that would occur if the boundary layer had been turbulent from 
the leading edge of the plate, and thus the bracketed term can be viewed 
as a correction that incorporates the effect of the actual location of 
transition on skin friction. Equation (30) is shown by figure 3. Note 
that a considerable increase in the friction factor may be expected if 
the transition Reynolds number is relatively high, but that the effect 
is small if transition occurs close to the leading edge of the plate. 


Equivalent Thermal Problem 

Reference 2 shows that the local Stanton number for an isothermal 
plate and a purely turbulent boundary layer is given by (neglecting 
temperature -dependent fluid-properties effects) 

St T (Re|jO)Pr 0 - 4 = 0.0296Re| 0 ' 2 (31) 

where 6 is measured from the virtual origin of the turbulent layer. 

The heat-transfer rate after an abrupt transition could be determined 
from this relation, with the Reynolds number based on the distance 
x-r, if it were not for the fact that the heat transfer over the lam- 
inar portion of the boundary layer differs from what would be trans- 
ferred to the turbulent layer in the distance from the virtual origin to 
the transition point. As a consequence, the condition of conservation 
of energy for the control volume aa’b’b (for dx small) cannot be sat- 
isfied if, in the equivalent thermal problem, the plate temperature is 
held constant. In order to allow the energy condition to be satisfied, 
it is necessary to adjust the plate temperature in the fictitious turbu- 
lent region before transition. It will be assumed that the plate temper- 
ature is constant over the fictitious turbulent region, and the value of 
this constant will be selected so that the total energy transferred to 
the laminar layer before transition equals that which would be trans- 
ferred to the fictitious part of the turbulent layer; that is, the part 
between the virtual origin and the real transition location. The heat 
transfer to the fictitious turbulent region may be determined from 


equation (31) , and the heat transfer to the laminar layer may be found 
from the -well-known Pohlhausen solution, 

St Zam Pr 2 / 3 = 0. 332Re” 0, 5 (32) 

The condition that the total heat transfers are equal may be written as 



AtGc p (o . 332Pr' 2 / ^e^ 0 ’ 5 jdx 



At'Gc. 


) (o.0296Pr' 0,4 Re^ 0,2 )di 


The indicated integrations lead to the result that 


At’ 

At 


1.012Pr 


-0.27 


(33) 


(34) 


Note that, for fluids with Prandtl numbers less than unity, the temper- 
ature difference in the hypothetical turbulent region between the vir- 
tual origin and the transition point is greater than the actual temper- 
ature difference, indicating that the total heat transfer in the laminar 
portion is greater than in the hypothetical turbulent region. This is 
possible because the fictitious turbulent region is much shorter than 
the actual laminar region. 


Local Stanton Numbers After Abrupt Transition 

The heat-transfer rates downstream of the abrupt transition may be 
determined by applying the methods of superposition described in refer- 
ence 4 to the equivalent thermal problem. Reference 4 shows that the 
local heat-transfer rate due to a series of steps in the wall temperature 
may be determined simply by summing the heat-transfer rates due to each 
step. Reference 3 shows that the heat transfer to a turbulent boundary 
layer for a plate that is unheated from the leading edge to l and main- 
tained at constant temperature thereafter is given by 

st_ [, m 9 / 10 

st T ' L “ 

where £ is measured from the leading edge. Thus, addition of the heat- 
transfer rate due to a step of height At' at the fictitious leading 
edge and the heat-transfer rate due to a step of height (At - At’) at a 
distance (s - r) downstream of the fictitious leading edge yields the 
following result: 
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St 


At = 0.0296Pr'°' 4 Re"°‘ 2 At 


x-r 


By using equation (28) to eliminate r , equation (36) may be reduced to 


+ 0.0296Pr'°' 4 Re"^ 2 jl - ( 7 * A 9 ^ 10 j (At - At') (36) 


St„ 


St < R V Ee s> r Ee s , s xT 0 ' 2 

(r^oT'L R%V s l 


X ll.012Pr' 0,27 + 


(l - 1.012Pr" 0 ’ 27 ^ 


X (1 - 


36.3 p. Be; 0 ' 375 
Re x s 


1 - (l - 36.3Bef- 375 ) 


9/10 


il/9 


= F(Re x ;Re s ) 


(37) 


where St(Re x ;Re g ) is the local Stanton number for the turbulent part of 

a boundary layer that is laminar from x = 0 to x = s and undergoes 
abrupt transition to a turbulent boundary layer at s; and Strp(Re x ;0) 

is the local Stanton number for a boundary layer that is turbulent from 
the leading edge, where x = 0 (as given by eq. (31) or an equivalent 
expression). Rote that, as Re g /Re x approaches zero, the local Stanton 

number approaches the value for a boundary layer that is entirely turbu- 
lent, and this is the correct limiting behavior. At the transition point 
the analysis indicates that the local Stanton number is and this is 

due to the discontinuity in the wall temperature in the equivalent thermal 
problem. Equation (37) is plotted for Pr = 0.7 in figure 4. Note that 
there is a considerable overshoot of the asymptotic curve (eq. (3l) ) if 
the transition occurs at a high Reynolds number; the local heat-transfer 
coefficient may be nearly twice as high as the values predicted by the 

D 

usual turbulent equation for Re g = 10 . At high transition Reynolds 

numbers the effect of transition location is very important and the in- 
crease in Stanton numbers persists well downstream of transition. It is 
probable that the increase in the local Stanton numbers and friction fac- 
tors is less if the transition occurs over a longer region, but it is 
evident that the effects of transition location are too large to be over- 
looked in many applications. The abrupt-transition analysis provides an 


4993 
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"upper limit" on these effects and is also useful in its own right in 
evaluating experiments where the boundary layer has been deliberately 
"tripped" in a two-d i mensional manner. 

RESULTS AND DISCUSSION 
Comparison of Abrupt-Transition Analysis 
with Heat-Transfer Experiments 

Using the flat-plate test equipment described in reference 2, a 
series of test runs were made in which abrupt transitions were stimulated 

5 5 

at Reynolds numbers ranging from 0.11x10 to 2.6X10 . In these runs 
transition was stimulated by a l/ 2-inch-wide strip of fine (414 grit) 
emery cloth, which was cemented to the plate at the desired point. A 
series of tests with different grits indicated that the 414-grit cloth 
had no appreciable "roughness effect" on the turbulent heat transfer 
downstream of the trip. In addition, one run was made in which the bound- 
ary layer was allowed to undergo natural transition. In this run transi- 
tion was not abrupt but started at a Reynolds number of about 2X10^ and 

appeared to be complete by a Reynolds number of 6X10^. The data from 
these tests are shown by figure 5 and are tabulated in table I. Refer- 
ence 2 shows that effects of temperature -dependent fluid properties may 
be lumped into a temperature ratio and that the local Stanton number for 
a purely turbulent flow may be determined from 

-0.4 

StrpPr 0,4 = 0.0296Re x °' 2 (38) 

where the fluid properties are to be evaluated at the free-stream temper- 
ature. Thus, the data are plotted in the form St(T w /T 00 ) < ^ > ’ 4 against 
Re x . Note that the data do tend to overshoot the limiting relation (38) 

and that the effect is greatest for the most delayed transitions. These 
are characteristics predicted by the foregoing abrupt-transition 
analysis . 

The experimental data from these runs are compared with the analysis 
in figure 6. The analysis appears to correlate the abrupt-transition data 
very well. The natural-transition data were corrected with the abrupt- 

transition analysis, based on a transition Reynolds number of 4X10 5 , which 
is approximately in the middle of the transition region. Only the fully 

turbulent data (Re x > 6X10 5 ) are shown in figure 6 for the natural- 

transition run. In spite of the fact that the transition is not abrupt. 
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the correction of the natural-transition data by the abrupt-transition 
analysis appears to be adequate. Thus, it seems that the abrupt- 
transition analysis may be used in estimating the effect of transition- 
point location for natural transitions if the Re s is evaluated some- 
where in the middle of the transition region. 


Comparison of Natural-Transition Heat-Transfer 


$ 

Data with Statistical Predictions 

c. 

By use of the abrupt- transit ion analysis (eq. (37)), the Pohlhausen 
laminar solution (eq. (32)), and the statistical prediction of the local 
Stanton numbers given by equations (ll) , the local Stanton numbers for 
the natural- transit ion run were predicted. The function 3(Re x ) was 

taken as an error integral, as is suggested by the results of Schubauer 
and Klebanoff (ref. 6) : 



Examination of the heat-transfer data for the natural-transition run in- 
dicates that transition begins to be important at about Re x = 2X10 5 and 
that the boundary layer is almost entirely turbulent by Re = 6X10^. 

Thus the mean, Re Sm , is about 4X10 . By taking Re^ as 10 D , the value 
of 3 is about 0.025 at Re x = 2X10^ and about 0.975 at Re x = 6X10^. 

As this is the approximate extent of the transition region. Re... was 

5 ^ 

taken as 10 in the predictions for the local Stanton number. The data 

are compared with the statistical predictions in figure 7. The predic- 
tions appear to be quite good, even in the transition region. Since the 
abrupt-transition analysis indicates that the local turbulent Stanton 
number at the transition is -oo, the predictions have singularities at 
Re 2<n ^ Res m ^ SSe eq - s ' ^ (37)). Since this is an 

artificiality introduced by the approximations of equations (10) , the 
predictions are not shown in the neighborhood of the singularities. The 
influence of the singularity extends only over a small region, however, 
as may be seen by examination of figure 4. Thus, it seems that the sta- 
tistical model, together with the abrupt -transition analysis and the dis- 
tribution function 3 , provides a method for calculating the heat trans- 
fer in the transition region and the effect of the location of transition 
on the turbulent heat-transfer rates if the location and approximate 
extent of the transition region are known. 


CONCLUDING REMARKS 


An analysis has been made that allows the prediction of the effect 
of the location of transition on the heat transfer to the turbulent bound- 
ary layer. The analysis assumes that the boundary layer undergoes an 
abrupt transition at some point downstream of the leading edge. The 
"patching” of the laminar and turbulent layers at the transition "point” 
is such that the momentum, continuity, and energy theorems are all satis- 
fied. The abrupt-transition data of the present investigation are in 
good agreement with the analysis. Unfortunately, these data were ob- 
tained over only a limited range of transition Reynolds numbers, and thus 
the analysis lacks experimental confirmation in the range where the ef- 
fect is largest. The abrupt-transition analysis may be used to predict 
the local turbulent Stanton numbers after a natural transition if the 
"transition Reynolds number" is evaluated in the middle of the transition 
region. If the extent of the transition region is known, it may be as- 
sumed that the function p(Re x ) is an error integral; then the statisti- 
cal treatment may be used to calculate the heat transfer to both the tran- 
sition and turbulent portions of the boundary layer. In this manner an 
estimate of the local heat-transfer coefficient can be made that is in 
good agreement with the experimental data for all types of flow. 


Stanford University, 

Stanford, Calif., October 22, 1957. 
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TABLE I 


EXPERIMENTAL HEAT- TRANSFER DATA 


G , 

lb 

°F 

1 w > 

Btu 

h r 

Btu 

St 

xio 3 

St © 0 ' 4 

* 10 3 

Re x 

XIO " 6 

G , 

lb 

<a 

q 

Btu 

V 

Btu 

St 

xlO 3 

fM 0 - 4 

st vd 

XlO 3 

Re x 

xlO - 6 

( hr ) ( sq ft ) 
* 10“ 3 

( hr)(sq ft ) 

( hr) (sq ft )(° F ) 

( hr) (sq ft ) 
« 10- 3 

( hr ) ( sq ft ) 

( hr) (sq f t ) ( ° F ) 

Natural 

transition ; t m = 

75 . 4 ° F ; p ,, = 0.0750 lb/cu ft 


Re 

a 

O . llxlO 5 ; = 65 . 3 ° P ; p ^ = 0.0765 lb/cu ft 


10.6 

23.9 

97.2 

4.07 

1.70 

1.73 

0.082 

~ 

11.7 

16.6 

165 

9.95 

3.53 

3.57 

0.092 

10.6 

24.4 

75.7 

3.10 

1.21 

1.24 

.136 

11.7 

16.3 

147 

9.04 

3.22 

3.26 

.151 

10.6 

24.6 

66.8 

2.72 

1.07 

1.08 

.186 

11.7 

16.3 

138 

8.48 

3.02 

3.06 

.206 

10.6 

24.2 

66.1 

2.73 

1.07 

1.09 

.236 

11.7 

16.3 

129 

7.94 

2.82 

2.88 

.264 

10.6 

24.7 

70.6 

2.86 

1.12 

1.14 

.286 

11.7 

16.2 

124 

7.65 

2.72 

2.76 

.319 

10.6 

24.8 

08.6 

3.56 

1.41 

1.44 

.334 

11.8 

16.1 

120 

7.48 

2.64 

2.68 

.377 

10.6 

24.5 

110.0 

4.49 

1.76 

1 . 79 

.385 

11.7 

16.3 

118 

7.24 

2.56 

2.59 

.433 

10.6 

23-8 

136.1 

5.77 

2.27 

2.31 

.434 

11.7 

16.3 

121 

7.40 

2.62 

2.67 

.489 

10.6 

24.2 

147.9 

6.11 

2.41 

2.45 

.485 

11.8 

16.4 

111 

6.74 

2.38 

2.41 

.546 

10.6 

24.1 

151.5 

6.29 

2.48 

2.52 

.533 

11.8 

16.4 

108 

6.61 

2.34 

2.37 

.601 

10.6 

24.0 

155.0 

6.46 

2.53 

2.58 

.587 

11.8 

16.4 

112 

6.86 

2.42 

2.45 

.660 

10.6 

23.9 

155.5 

6.51 

2.55 

2.60 

.637 

11.8 

16.7 

110 

6.57 

2.32 

2.35 

.715 

10.6 

24.2 

156.7 

6.48 

2.55 

2.60 

.684 

11.8 

16.6 

107 

6.46 

2.23 

2.31 

.771 

10.6 

24.2 

151.5 

6.26 

2.61 

2.66 

. 732 

11.8 

16.4 

109 

6.67 

2.35 

2.38 

.827 

10.6 

24.4 

156.4 

6.37 

2.51 

2.55 

.784 

11.8 

16.4 

107 

6.54 

2.31 

2.34 

.884 

10.6 

24-7 

150.9 

6.11 

2.55 

2.59 

.831 

11.8 

16.9 

108 

6.38 

2.26 

2.29 

.939 

10.6 

24.5 

144.3 

5.89 

2.32 

2.36 

.883 

11.8 

16.9 

105 

6.20 

2.19 

2.22 

.995 

10.6 

24.4 

139.6 

5.78 

2.28 

2.32 

.933 

11.7 

16.9 

102 

6.01 

2.13 

2.16 

1.048 

10.5 

24.6 

141.5 

5.75 

2.27 

2.32 

.976 

11.8 

16-7 

101 

6.07 

2.14 

2.17 

1.104 

10.5 

25.0 

141.8 

5.67 

2.24 

2.29 

1.029 

11.7 

16-8 

104 

6.17 

2.19 

2.22 

1.160 

10.5 

25.3 

135.5 

5.36 

2.12 

2.16 

1.078 

11.7 

16.6 

93 

5.58 

1.98 

2.02 

1.214 

10.5 

25.1 

132.6 

5.29 

2.09 

2.13 

1.127 

11.7 

16.9 

99 

5.84 

2.07 

2.10 

1.270 

Re s 

= 0 . 92 * 10 5 ; t m = 

69 . 9 ° F ; p m - 0.0755 lb/cu ft 


Re s 

1 . 44 x 10 5 ; t . = 71 . 3 ° F ; p „ - 0.0756 lb/cu ft 


11.8 

16.3 

96.5 

5.92 

2.10 

2.12 

0.092 

11.7 

16.7 

72.7 

4.35 

1.55 

1.56 

0 . 092 

11.8 

16.1 

171.3 

10.64 

3.77 

3.82 

.151 

11.7 

17.5 

50.9 

2.91 

1.04 

1.05 

.150 

11.7 

16.1 

137.7 

8.55 

3.04 

3.08 

.206 

11.7 

16.3 

159.2 

9-77 

3.48 

3.53 

.205 

11.8 

16.1 

131.8 

8.18 

2.88 

2.92 

.264 

11.7 

16.5 

149.2 

9.04 

3.22 

2.36 

.261 

11.8 

16.2 

125.4 

7.75 

2.73 

2.78 

.320 

11.7 

16.7 

138.9 

8.32 

2.97 

3.00 

.316 

11.8 

16.1 

120.2 

7.48 

2.64 

2.67 

.376 

11.7 

16.7 

132.7 

7.95 

2:82 

2.86 

.373 

11.8 

16.2 

116.9 

7.21 

2.56 

2.59 

.429 

11.7 

16.6 

124.6 

7.50 

2.67 

2.70 

.427 

11.7 

16.3 

121.0 

7.42 

2.64 

2.67 

.485 

11.7 

16.6 

127.9 

7,70 

2.74 

2.77 

.484 

11.8 

16.4 

112.0 

6.83 

2.42 

2.45 

.542 

11.7 

17.0 

117.1 

6.89 

2.45 

2.48 

.540 

11.8 

16.4 

109.8 

6.70 

2.36 

2.39 

.600 

11.7 

16.8 

116.0 

6.90 

2.45 

2.48 

.595 

11.7 

16.5 

113.1 

6.87 

2.44 

2.47 

.652 

11.7 

17.0 

115.8 

6.82 

2.42 

2.45 

.650 

11.8 

16.7 

111.1 

6.66 

2.35 

2.38 

.712 

11.8 

17.2 

117.4 

6.83 

2.41 

2.44 

.710 

11.8 

16.4 

107.5 

6.56 

2.32 

2.36 

.765 

11.8 

17.2 

113.8 

6.62 

2.35 

2.38 

.765 

11.8 

16.5 

103.6 

6.27 

2.21 

2.24 

.824 

11.7 

17.5 

109.6 

6.25 

2.21 

2.24 

.820 

11.8 

16.5 

103.9 

6.30 

2.22 

2.45 

.881 

11.8 

17.1 

112.6 

6.58 

2.23 

2.35 

.880 

11.8 

16-5 

105.0 

6.36 

2.24 

2.27 

.938 

11.8 

17.3 

109.3 

6.32 

2.23 

2.26 

.935 

11.8 

16.5 

102.7 

6 , 22 

2.19 

2.22 

. 392 

11.6 

I t .3 

111. 5 

6.44 

2.27 

2.30 

.991 

11.8 

16.6 

99.7 

6.01 

2.13 

2.16 

1.042 

11 .8 

17.4 

108.1 

6.21 

2.19 

2.22 

1.048 

11.7 

16.6 

102.3 

6.17 

2.19 

2.22 

1.092 

11.7 

17.4 

104.4 

6.00 

2.13 

2.16 

1.093 

11.8 

16.8 

102.0 

6.08 

2.15 

2.18 

1.156 

11.7 

17.5 

107.6 

6.14 

2.18 

2.21 

1.150 

11.8 

16.7 

94.3 

5.65 

2.00 

2.03 

1.210 

11.7 

17.7 

99.8 

5,61 

1.99 

2.01 

1.208 

11.8 

16.0 

93.7 

5.58 

1.98 

2.00 

1.266 

11.7 

18,0 

103.0 

5.77 

2.05 

2.07 

1.264 

Re g - 2 . 00 * 10 5 ; t _ - 71 . 5 ° P ; p m « 0.0758 lb/cu ft 

Re = 
s 

2 . 60 * 10 5 ; t . » 64 . 0 ° F ; p . = 0.0766 lb/cu ft 


11.7 

17.8 

76.2 

4. 38 

1.53 

1.55 

0.091 

11.7 

18.2 

78.0 

4.28 

1.52 

1.55 

0.092 

11.8 

18.4 

58.2 

3.16 

1.12 

1.14 

.150 

11.7 

18.3 

60.9 

3.33 

1.17 

1.21 

.151 

11.8 

10.6 

38.7 

2.08 

.74 

.75 

.206 

11.7 

18.1 

53.6 

2.96 

1.05 

1.07 

.208 

11.8 

17.3 

165.4 

9.57 

3.39 

3.44 

.262 

11.7 

18.3 

50.5 

2.76 

.98 

.99 

.264 

11.8 

17.4 

151.0 

8.68 

3.07 

3.16 

.317 

11.7 

17.9 

153.0 

8,55 

3.04 

3.09 

.316 

11.8 

17.7 

144.0 

8.13 

2.87 

2.91 

.373 

11.7 

17.9 

155.2 

8.68 

3.08 

3.12 

.375 

11.8 

17.8 

138.9 

7.80 

2.76 

2.80 

,428 

11.7 

18.1 

147.5 

8.15 

2.89 

2.94 

.432 

11.0 

18.1 

135.5 

7-48 

2.55 

2.56 

.484 

11.7 

18.1 

146.3 

8.09 

2.88 

2.92 

.487 

11.9 

18.1 

122.2 

6.76 

2.38 

2.41 

.544 

11.7 

18.5 

134.5 

7.27 

2.58 

2.62 

.544 

11.8 

17.9 

120.7 

6.74 

2.37 

2.41 

.598 

11.7 

18.3 

130.7 

7.14 

2.54 

2.52 

.600 

11.8 

18.0 

121.2 

6.74 

2.36 

2.41 

.656 

11.8 

18.5 

131.3 

7.10 

2.51 

2.55 

.662 

11.8 

18.3 

120.4 

6.58 

2.32 

2.36 

.711 

11.7 

18.6 

130.2 

7.00 

2.49 

2.53 

.710 

11.8 

16.1 

118.4 

6.55 

2.32 

2.35 

.764 

11.7 

18.6 

129.5 

6.96 

2.47 

2.51 

.768 

11.8 

18.3 

116.4 

6.36 

2.24 

2.27 

.823 

11.7 

16.5 

123.0 

6.65 

2.36 

2.40 

.824 

11.9 

18.0 

111.8 

S .21 

2.18 

2.22 

.881 

11.7 

18.4 

122.8 

6.67 

2.37 

2.41 

.880 

11.8 

18.2 

113.2 

6.23 

2.19 

2.22 

.935 

11.7 

18.4 

118.9 

6.46 

2.30 

2.33 

.935 

11.8 

18.4 

110.6 

6.01 

2.13 

2.16 

.985 

11.7 

18.5 

116.2 

6.28 

2.13 

2.16 

.990 

11.8 

10.2 

106.0 

5.83 

2.06 

2.09 

1.040 

11.7 

18.4 

112.6 

6.11 

2.17 

2.20 

1.049 

11.8 

18.1 

103.7 

5.72 

2.03 

2.06 

1.092 

11.7 

18.4 

109.7 

5.95 

2.12 

2.14 

1.110 

11.8 

18.2 

105.7 

5.80 

2.06 

2.08 

1.151 

11.7 

18.3 

112.5 

6.15 

2.18 

2.20 

1.160 

11.8 

18.4 

102.1 

5.55 

1.96 

1.99 

1.210 

11.7 

18.2 

102.6 

5.64 

2.00 

2.03 

1.216 

11.8 

18.7 

104.8 

5.60 

1.98 

2.01 

1.262 

11.7 

18.5 

105.3 

5.70 

2.01 

2.06 

1.270 





Figure 1. - Transition model 







Figure 3. - Effect of delayed transition on turbulent friction factors. 







Figure 5. - Effect of location of transition on turbulent heat transfer. Prandtl number, 0.7 
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Figure 6. - Comparison of data with abrupt- 



analys i s . Prandtl number , 0.7. 
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Figure 7. - Comparison of natural- transition data with statistical predictions. Prandtl 
number, 0.7. 
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